Boundary Term Contribution to the Volume of a 
Small Causal Diamond 



Surbhi KhetrapaP and Sumati Surya b , 
a BITS Pilani K. K. Birla Goa Campus 
b Raman Research Institute, Bangalore, India 

December 5, 2012 



Abstract 

; 

^-l In his calculation of the spacetime volume of a small Alexandrov inter- 

1 ■ val in 4 dimensions Myrheim introduced a term which he referred to as a 
£^ surface integral [1] . The evaluation of this term has remained opaque and 

led subsequent authors to evaluate the volume using other techniques [2] . 
It is the purpose of this work to demystify this integral. We point out that 
it arises from the difference in the flat spacetime volumes of the curved 
and flat spacetime intervals. An explicit evaluation using first order de- 
generate perturbation theory shows that it adds a dimension independent 
factor to the volume of the flat spacetime interval as the lowest order cor- 
rection. Our analysis admits a simple extension to a more general class 
of integrals over the same domain. Using a combination of techniques we 
also find that the next order correction to the volume vanishes. 

1 Introduction 

Spacetime volume plays a crucial role in the discrete-continuum correspondence 
of causal set quantum gravity [3J. In this approach the fundamental entity 
underlying the continuum is a locally finite partially ordered set or causal set. 
The order relation of the causal set corresponds to the causal structure of the 
spacetime, while the local finiteness means that underlying every A" Planck 
volumes of a spacetime region there are, on average, A" of elements of the causal 
set. Calculations of the volumes of generic but causally well-defined spacetime 
regions are therefore of particular interest to causal set theory. 

In an insightful CERN preprint on discrete statistical geometry Myrheim [T] 
introduced several key concepts that were subsequently adopted by the causal 
set approach [3J. An important calculation in this work is of the volume V of 
a small Alexandrov interval I[p, q] between two chronologically related events 
p and q. The metric is expanded about the mid-point r of the flat spacetime 
geodesic 70 from p to q using Riemann Normal Coordinates (RNC), with the 
smallncss parameter being the proper time T from p to q. In 4 spacetime 
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dimensions Myrheim obtained an expression for V to order T 2 which depends 
on the scalar curvature as well as the curvature components along 70. While 
the calculation is straightforward to set up, the evaluation of one of the terms, 
referred to as a "surface integral" in pQ has remained opaque. This led Gibbons 
and Solodukhin to an evaluation of V which side steps this integral by instead 
calculating the universal coefficients that appear in an order T 2 expansion of 
the volume for the Einstein Static Universe and de Sitter spacctime [2]. 




P 

Figure 1: An Alexandrov interval I[p, q] in flat spacetime. T is the proper time 
between the events p and q. 

It is the main purpose of this work to decode Myrheim's "surface" integral. 
Rather than viewing it as a surface term, we find it more fruitful to think of 
it as the difference in the flat spacetime volumes of the regions defined by the 
curved spacetime interval I[p, q] and the flat spacetime interval Io[p, q\- As we 
will show in Section 2, the boundary of I[p, q] can be determined to order T 2 
by merely considering the modification to the light cones in the tangent spaces 
T p M and T q M of p and q respectively; the effect of the acceleration of the 
null geodesies can be shown to be sub-leading. The modification in T p M can 
be understood from simple first order degenerate perturbation theory. In the 
RNC this curvature dependent perturbation is in fact restricted to the spatial 
directions, thus generically transforming the spherical cross section of the light 
cone to an ellipsoidal one. Making a coordinate transformation to the principal 
directions of this ellipsoid we find an expression for the sum of the modified 
principal eigenvalues using first order degenerate perturbation theory. We show 
that this sum appears crucially in the Jacobian of transformation relating the 
flat spacetime integral to Myrheim's term to lowest order thus completing the 
first part of our analysis. In Section 3 we evaluate the Myrheim term and show 
that it provides a dimension independent factor to the flat spacetime volume to 
this order. This allows us to complete the RNC calculation of the volume V of 
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I\p, q] up to order T 2 for arbitrary n. Our expression agrees with that obtained 
by Myrheim [1] for n = 4 as well as that obtained by Gibbons and Solodukhin 
[2] for general n using a different technique. An important by-product of our 
analysis is that it allows us to calculate more generic integrals which crop up in 
causal set theory. We give an example of one such integral which has recently 
been used to obtain a new expression for the scalar curvature of a small causal 
set [5] using a curved spacetime generalisation of the work of [6]. 

What of higher order corrections to the volume? In [7] this question was 
partially addressed in the discussion on the volumes of large causal diamonds 
and their asymptotic behaviour. As evident from the RNC expansion, higher 
order corrections to the volume will include higher derivatives of the curvature. 
In Section 4 we will examine these terms using the techniques developed above. 
While the light cones in T p M and T q M can be studied using the next order per- 
turbation analysis, the acceleration of the null geodesies is non-negligible. Hence 
it does not suffice to look at T P M and T g M, although it does help determine the 
types of higher derivative terms appearing to this order . Using a combination 
of the analysis of light cones in T p M, T q M and the approach of Gibbons and 
Solodukhin [2] we evaluate the volume to 0(T 3 ) in FRW spacctimcs. We find 
that the 0(T 3 ) correction to the volume vanishes altogether. 

Because of the appearance of higher derivatives of the curvature it is tempt- 
ing to ask if the volume expansion can have any significance in determining 
effective actions in quantum gravity. This idea is not so far-fetched in causal set 
theory where spacetime volume plays a fundamental role, making it possible to 
speculate that higher derivative corrections to the action must be determined by 
these corrections to the volume. This could distinguish the causal set approach 
from other approaches to quantum gravity. While our primary focus in this 
note is a modest one, namely to decode Myrheim's calculations and to seek the 
extension of that analysis, it can be seen as a first step towards this more ambi- 
tious goal. The vanishing of the 0(T 3 ) term then suggests that the corrections 
to the Einstein-Hilbert action in the effective continuum action arising from 
causal set theory will only arise from second derivatives of the curvature, mod- 
ulo boundary terms. However, a systematic approach to prove such a conjecture 
is currently beyond our scope. 

2 The Boundary of I[p,q] 

The RNC about a point r £ M is defined within a convex normal neighbourhood 
Q of r in the spacetime (M,g), i.e., a region Q C M in which the exponential 
map exp : T p M — > Q is a diffeomorphism for any p S Q. In RNC, the geodesies 
emanating from r are used to coordinatise Q, and the spacetime metric at r is 
taken to be flat, i.e., <? a b(i") = Vab- This implies that the Christoffcl connection 
rjj b (r) = 0, so that the metric at any x € Q can be expanded as 

g ab {x) = i M {0) - ^x c x d R acbd {0) + 0(x 3 ), (1) 
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where we have used the RNC identity d^tc) (0) = 0. Here, and in the future 
we will often use the short form to denote the origin r = (0,0, ...0) of the 
RNC. 

The volume of a small Alcxandrov interval I[p, q] in n spacctime dimensions 
between the points p = (-T/2, 0, 0, . . . , 0) and q = (T/2, 0, 0, ... 0) in RNC is 
therefore given by the integral 

V = JV=g d n x = J (l - ^x c x d R cd (0) + 0(x 3 ^j d n x. (2) 
ib,<?] ib,?] 

which was first calculated for n = 4 by Myrheim [1] . Importantly, the integration 
is over a region I[p, q] which itself is determined by g a b{x) and hence contains 
corrections to the flat spacetime interval I [p, q]. Up to 0(x 2 ) this integral can 
be split into two parts V = 2m + 2 2 , where 

Im=Jd n x, I 2 = J (-^x c x d R cd (0)^jd n x. (3) 
i[p.?] iob,g] 

As we will show in the next section, the second term is straightforward to 
evaluate in arbitrary dimensions, much of the simplification arising from the 
fact that odd terms do not contribute to lo\p, q] because of its symmetries. 

The integral 2m itself comprises two pieces: the volume Vo of Io[p, q] plus 
a contribution lAi[p.q] which Myrheim referred to as a "boundary" term. This 
term was evaluated for n = 4 in [T] without any explanation, although how to 
evaluate it it has been far from obvious to subsequent researchers [2J. Indeed, it 
is this precise integral that we wish to decode in the present work. We first note 
that this term is not a boundary integral but simply the difference in the flat 
spacetime volumes of the interval I[p, q] with respect to rj a b and g a b respectively, 
as realised also in [2j. 

Thus, in order to determine 2m, we need to first obtain the boundary of I[p, q] 
since the difference in the volumes lAi[p,q] roughly arises from the difference in 
the range of integration. In flat spacetime the boundary of Io [p, q] is that of a 
pair of uniform (light)concs with base radius T/2, one emanating to the past 
from q and the other to the future from p so that we can evaluate 

T/2 T/2-t 

Vo = / d"x = 2 [dt [ r- a dr / dfi„_ 2 = (4) 



I [p,q] S' 
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where A„_2 is the volume of a uniform n — 2 sphere S n ~ 2 . 

The presence of curvature obviously modifies the boundaries of the two light 
cones, and it is this that we now attempt to quantify. The boundary of I[p, q] is 
ruled by future-directed and past-directed null geodesies emanating from p and 
q respectively. The tangents to these null geodesies in turn lie along the future 
and past light cones in T p M and T q M respectively. The effect of curvature in 
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general would be to accelerate these null geodesies, so that they no longer lie 
along T p M once they leave p, and similarly for q. Thus, one would no longer 
expect as simple an integration as in Eqn. 

Nevertheless, the first step to take is to determine the future and past light 
cones in T p M and T q M , respectively Because of the symmetry of Io[p, q] it 
suffices to restrict our attention to T q M. Using the RNC expansion of the 
metric at q 

gab(q)=Vab(0)-^T 2 R Qa0b (0), (5) 

the tangents £ a = = (C°, C 1 ! • ■ • C" 1 ) to the past and future directed null 
geodesies at q satisfy 

9ab(q)CC b = 



-(C ) 2 + £(O 2 = ^T 2 iW0)CV (6) 

i=l 

where we have used the symmetries of the Riemann tensor to simplify the ex- 
pression. Thus the light cone in the tangent space T q M is given by the matrix 
equation 

C*MC = (C°) 2 , (7) 
where £ is the spatial part of ( a and 

= Sij - T 2 —Rom 3 - ( 8 ) 

When there is no curvature Eq (JT]) reduces to (£) 2 = (C ) 2 which is the equation 
for a uniform (light)cone, i.e., a cone with (n — 2)-dimensional spherical cross- 
sections with radii C°- 

Since M is a real symmetric matrix, it can be diagonaliscd by an orthogonal 
matrix Q, so that inserting Q T Q = 1 in Eq yields 

fA|*=(C ) 2 =* ^A,(e) 2 = (C ) 2 , (9) 

i 

where ^ = QC and A is the diagonalised form of M, with Ay = 5yA, This 
gives an equation for a spatial ellipsoid for fixed £° with principal axes along 
the components of £. Thus the presence of curvature modifies the regular light 
cone in T q M into one whose spatial sections are ellipsoids with principal axes 
given by the £ l . 

A past-directed null geodesic from q with tangent vector £ a = (C°,£) satis- 
fying Eq © obeys the geodesic equation 

^i* = -n c m b u% = -!« c (o)eV c i 9 , (io) 

in RNC. In other words, to lowest order, the acceleration of geodesies is zero, 
as expected. Expanding along the affine parameter s 

C = C\ q s^d TU0)e\ q t% + 0(s 2 ), (11) 
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we see that for small s, the second term is sub- leading to order 0(T 2 ) but not 
to order 0(T 3 ). Thus, if we find that the first correction to £ a is of order 0(T 2 ), 
then the acceleration is sub-leading, so that the boundary of I[p, q] is determined 
by the light cones in T PjQ M alone, to order T' 2 . 

Wc find the leading order correction to £ a by parallel transporting it from 
r = (0, 0, . . . , 0) to q along a geodesic. Indeed, as is easily shown, up to order 
T 2 the curve 7 = (t, 0, . . . , 0) from r to q is a geodesic. The tangent to 7 is 
U a = (1, 0, . . . , 0), so that dolI b = all along 7. For 7 to be a geodesic therefore 

U a V a U\t) = d U b (t) + T b Q (t) = td o T b O0 (0) + 0(t 2 ) (12) 

must vanish to the required order. 
Using the RNC identity 

d a lt(0) = -r s (R abc d (0) + R a J(0)) (13) 

and the properties of the curvature tensor, we see that U a V a U b (t) = 0(t 2 ), so 
that 7 is a geodesic up to 0(T 2 ). 

Now consider the parallel transport of a vector us a from r to q along 7, 
U b V b uj a = : 

d Q uj a + td r a ob (o)u b = => d oj a - ^i? ofco Q (o) w " = (14) 

where we have used the identity Eq (|13|) . The symmetries of the Riemann 
tensor imply that the time-component of any vector does not change along 7, 
i.e., <9 £° = up to this order. Since i? jO (0) ls symmetric and real (and hence 
Hermitian) let us use its eigenfunctions {e^ 1 } R o ^(0)e^ k ^ = p^e^ 1 to give 
us a (spatial) orthonormal basis vectors. Thus the parallel transport equation 
reduces to 

d e^ l R ^(Q)e^ = 
e^ z (t) =e (fc)? (0)cxp^"c=) 

e^CT/2) « U + ^ V(0))e (fc)i (0) (15) 

to this order. In other words, the spatial components of any vector have an order 
0(T 2 ) correction from fiat spacetime. Thus, we have shown that the acceleration 
of the null geodesies at q is negligible to this order and can therefore be ignored. 

We thus find that to 0{T 2 ) the modified past light cone from q is identical to 
that in T q M and hence, generically has ellipsoidal instead of spherical sections 
at constant time. In order to evaluate the boundary of q] it therefore suffices 
to find the relevant properties of this light cone in T q M, in particular to find 
the eigenvalues Xi of M . 

The Xi can be determined to order T 2 using first order perturbation theory 
in T 2 , or equivalently, second order perturbation in T where the first order 
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perturbation is zero. We prefer to use the latter terminology since it will make 
it easier to adapt to the 0(T 3 ) perturbations. Rewriting 

M = M {n) + T 2 M {2) (16) 

where M^> = I and M-^ = — j^Rowj, we note that the zeroth order eigenfunc- 
tions are simply the unit vectors i/'S in the principal directions, with eigenvalues 
\[°) = 1. The first correction A. 2 ? to the x[°) are then the eigenvalues of the 

operator ij;^ M^ip^, which in this case is simply itself. 

It is in fact not necessary for us to solve the exact eigenvalue problem in order 
to obtain an expression for the volume up to this order. As will become apparent 
shortly, all that is required is to obtain the sum A^ ; . The characteristic 
equation ||M( 2 )-A/|| =0 gives rise to an (n — l)th order polynomial 

A n-1 _ yv-2 A4 2) + 0(X n ~ 3 ) = 0. (17) 

(2) 

Rewriting this equation in terms of its roots A^ 

n-l , s 
i ^ ' 

we see that the coefficient of the (n — 2)th order term is 

-E A S=-E M i? )= E = (19) 

i i i 

(2) 

While we will not require for our main purpose to solve for X)J, it is never- 
theless instructive to see the effect of curvature on the light cones in a simple 
example. For n = 3 

AL 2) = \ (Mi (2) + M$ ± y/(M™-M™ ) 2 +4(A4 2) )^ . (20) 

Thus, the light cone in T q M is transformed from a symmetric cone to one with 
elliptical sections at constant £° with the lengths of the semi-major and semi- 
minor axes being £° / yfxZ and C°/y / A+, respectively. The principal directions 

(2) (2) (2) (2) (2) 

are degenerate, i.e., A+ = A_ iff M^' = M 22 ' and M[ 2 ' = 0, i.e., the curvature 
is isotropic. In this case A± = 1 — ^jT 2 Rqo, which means that the light cone, 
though still symmetric, is narrower than in fiat spacetime if i?oo > and wider 
otherwise. We can expect a straightforward generalisation of this behaviour in 
higher dimensions - the light cone could be wider in some spatial directions 
and narrower in others depending on the curvature. These various possibilities 
are already evident in n = 3. For example, if Yli-^u = — j^Rqo is equal to 
A = ((M} 2) - M 2 (2) ) 2 + 4(M 1 (2) ) 2 )V2 (i.e., M 12 = y/M u M 22 ), A+ < 1 for 
Roo > and A + > 1 for R no < while A_ = 1 in both cases. If i?oo > 
and A > j^R QOl then A+ > 1 and A_ < 1, etc. Figure [2] shows one of these 
possibilities. 
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Figure 2: The effect of curvature on the boundary of I[p,g] up to order T 2 . The 
dotted lines indicate the symmetric boundary with respect to flat spacetime and the 
bold lines the modified boundary in the presence of curvature. In this figure the two 
principal directions are taken to be non-degenerate. 

3 The Volume Calculation 

We are now in a position to calculate the volume of I[p, q] in arbitrary dimen- 
sions. But first let us compare this approach with the one in [2J. The starting 
point in [2J was to assume a form for the volume 

V = V (l + a(n)R(0)T 2 + p{n)R QQ {0)T 2 + 0(T 3 )), (21) 

based on Myrheim's calculation, thus bypassing the explicit evaluation of Eq 
l[2"]). The universal constant a(n) was evaluated by calculating the volume of a 
causal diamond explicitly for the Einstein static universe for which Rqq is iden- 
tically zero, and expanding in powers of T . Using this, the universal constant 
P{n) was calculated from the volume of a causal diamond de Sitter spacetime 
(R,Roo 7^ 0), again by expanding in powers of T. Thus, the authors of [2 
obtained a general formula for the volume of a small causal diamond in n space- 
time dimensions to order T 2 , and this will serve as a check for our calculations. 
Moreover, this technique of [5] will prove to be useful in our exploration of the 
next higher order corrections to the volume in Section 4. 

As shown in the previous section, the contribution from /m to the volume 
to 0(T 2 ) comes from a region bounded by light-cones whose tangent vectors at 
q satisfy Eq (|9]). The boundary of the backward light-cone from q is therefore 
given by 

(^-'l^EW, (22) 

i 

where we have chosen the spatial coordinates j/j along the principal axes of the 
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spatial ellipsoids. We may therefore construct the nested integral 

Im = 2 / dz / dy 1 dy 2 --- / dy n -\ (23) 



— wo — Wi 



u (»-2) 



where zq = T/2 — t, Xk+iWk = y z 2 — Xh=i \i)Vi- The coordinate trans- 
formation Z{ = \f\iyi simplifies the integral to that of a regular (uniform) cone 
along with a factor corresponding to the Jacobian of transformation 

1 = n x x = 1 T E A « + ° (r3) ' (24) 

since A t = 1 + T 2 A[^ + 0(T 3 ). Thus, 

= V (l — J2 A (i)J = Vo (^1 + ^oo(0)J (25) 

where we have used Eq (|19p . Thus, our main result is that in any dimension 
the boundary term of Myrheim [T] takes the form 



/aim = V I — i2oo(0) J. (26) 

Wc complete this section by including the computation of the I 2 integral in 
arbitrary dimensions. Let us rewrite 



io[p,?] 

where 



h= J (-^x c x d R cd (0)jd n x = -^(J 1 + J 2 + J 3 ) 



(27) 



/n— 1 ,, n— 1 ,, 

tTa: t 2 i?oo(0), J 2 = 2 ^ / t^x t i? *(0), J 3 = J] / d"x xV i? y 

Io[p,«] l_1 Io[p,9] Iobi9] 

(28) 

Evaluating 



Ji = 2^00(0) J dtt 2 J dr r n ~ 2 J dQ n 
s n 

T 
2 

i? 00 (0) Jdtt 2 (^--t 



g „-a 
ra-1 



n- 1 



n(n-l)(n+l)(n + 2) 





n+2 

«oo(0)(^) • (29) 
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To evaluate J2, J3 it is convenient to transform from the Cartesian coordinates 
x % to spherical polar coordinates x l = rf l (£l), where 



n, 02 ■ ■ ■ V n -2 



) = 



n—i—1 

1] sin0 fe cos0 n _i (i > 1) 

k=l 

n-2 

rising (i = l) 

fe=l 

Since the angular contribution to J2 is odd, it must therefore vanish. Similarly, 
the only contribution to J3 comes from the even terms. The contribution from 
each x l being the same, by spatial spherical symmetry of Io[p, q] , we may simply 
evaluate the integral for x n ~ l = rcos0i 



n—l ~ 

J3 = ]T / d n x (x*) 2 R u (0) 



lob, 1 
n-X 



— 4--t 



= 2^ J R ii (0)^ J dt J drr n J dfl n - 2 cos 2 
1=1 S"- 2 



(30) 



where 



J dfi„_ 2 cos 2 0i = j <Wi...<W n _a^JJ sin("- 1 - 2 ) 0^ cos 2 0i 



so that 



Thus 



d6»i sin™" 3 6»i cos^ 9 1 I d9 2 sin"" 4 



n— 4 , 



r/0 



n-2 



J 3 = 





n — l 



2A„_ 



(n-l)(n + l)(n + 2) 



n—l \ /rp\ n+2 



(31) 
(32) 



1 



3(n- l)(n + l)(n + 2) V 2 
nV 



Tl + 2 



- + 1 i? oo (0) + i?(0) 



24(n+l)(n + 2) \\n 
Adding this to Immi we get 



T - + 1 i? O o(0) + i?(0) 



(33) 



V = V 1 



24(n + 1) 



T 2 i?oo(0) 



24(n + l)(n + 2) 



T 2 i?(0) (34) 
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which matches Myrheim's expression for n = 4 [T] as well as the expression in 
[2] for arbitrary n. 

As discussed in the introduction, an important reason to be able to calcu- 
late the integral 7m explicitly is so that we can extend the analysis to other 
integrals that make their appearance in causal set theory. In [5], for example, 
the following integral must be evaluated for integer m > 0: 



K= I d n y T "\ (35) 

I[p>9] 

where r is the proper time from the event that one is integrating over and the 
future most point q of I[p, q]. Since the boundary of the light cones are given 
by Eq @, we see that r 2 = (T/2 — t) 2 — ^ Aij/f . As in the evaluation of 7m 
the limits of the integrals are again given as in Eq (|23[) . Again, performing a 
change of coordinates = s/Xlyt simplifies the integral to 

K = — 1= J d n z r™ = (l+ ^i?oo(0)) J d n z t™. (36) 
lob,?] io[p,<?] 

The universality of this term proves to be crucial in determining a recurrence 
formula for the average numbers of "fc-chains" in a causal set C which is ap- 
proximated by q]. This in turn allows us to find an expression for the scalar 
curvature as well as the dimension in purely order theoretic terms [5]. 

4 On the Order T 3 Corrections to the Volume 

In this section we show that the next higher order correction 0(T 3 ) to the 
volume is in fact zero. 

To order T 3 , the volume is given by 

V= J "V=fl d n x =J(l - ^x c x d R cd (0) - ^x e x c x d d e R cd (0) + 0{x i )^jd n x. 
i[p>«] i[p>?] 

(37) 

where we have used the 0(x 3 ) expansion in the RNC 

9ab(x) = Vab (0) - lx c x d R acbd (0) - \x e x c x d d e R acbd {Q) + 0{x 4 ) (38) 
3 6 

which arises from the identity d^ a dbT^(0) — or equivalently 

d d d {e T a b)c (Q) = ~ \2d {e R a b)cd (0) + d e R a cbd (0) + d d R a ecb (Q)^j . (39) 

Again, we may split up the above integral into a piece that comes entirely 
from the flat spacctime interval Io[p, q] and a "remainder" coming from the 
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difference in the two integrals 

V = h+h + h 
h = J d n x 
i[p,s] 

^x c x d R cd (0)d n x, 
ib,?] 

h = j (-^x e x c x d d e R cd + 0(x 4 )Jd n x. (40) 
lob,?] 

The symmetry of Io [p, q) means that the leading order contribution to Z3 , being 
odd, vanishes. In addition, since the leading order correction to I[p, q] is of 
order T 2 and the integrand of I2 is 0(x 2 ), this term must be 0(T 4 ). Hence the 
only 0(x 3 ) contribution can come from l\ which is analogous to the Myrheim 
integral. 

The 0(T 3 ) corrections to V should be of the form V a W b Z c d a R bc (0) where 
the V a , W a , Z c can be time-like or space-like vectors. Since to 0(T 2 ) one only 
has terms of R and i? o, a first guess is to include terms like d R and <9 i?oo- 
We show that these are precisely the sorts of terms that appear when examining 
the light-cones in T p M and T q M to this order. 

To order T 3 , the metric at q is 

g ab (q) = Vab(0) - -^T 2 R 0a0b (0) - ^T 3 d R aOb{0), (41) 
so that for £ a G T q M satisfying the null condition g a b(q)( a ( b = 

- (C°) 2 + E(C) 2 = (^T 2 Ro i0j (0) + ±T 3 d R 0i0j (0))et j , (42) 



1=1 



or ( M( = (C ) 2 which is the next order modification of Eq (J7J), where now 
Mij = 5 tj - ^T 2 Ro i0j - ^T 3 d Ro i0j = + T 2 M^ + T 3 M^ . (43) 

This suggests that only the time-derivatives of the Riemann tensor is relevant 
to the volume calculation. Moreover, using Eqn. (|39|) we see that the curve 7 
from r to q along the i-axis is also a geodesic to 0(T 3 ). This follows from the 
fact that 

U a V a U b = jd 2 TU0) + 0(t 3 ) = 0(t 3 ), (44) 

since 9grg Q (0) = from ((39)) . Moreover, again to this order the time-component 
of any vector remains unchanged under parallel transport while the spatial com- 
ponents satisfy 

+2 



d e^ - ^iV(0)e<^ - L do R OjO \0)e^ = 0. (45) 
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Thus, at least in T q M, even to order T 3 , the only derivative of the curvature 
that appears is the time derivative. If we were to ignore the acceleration of the 
null geodesies, and use third order perturbation theory to find the eigenvalues 
of M, then again, I\ differs from the flat space-integral purely by the Jacobian 
J Eq pif . which to this order is 

J=l-^-£Ag-T3£Af. (46) 

i=l i=l 

n-1 

Given the form of M^ 3 \ this suggests that only the trace ^2 doRoMj = 9qRoo 

i=i 

will contribute to J. We show that this is indeed the case, using a simple exten- 
sion of standard third order perturbation theory. Expanding M, its eigenvalues 
and its eigenfunctions in powers of T 

(0) + TM {1) + T 2 M {2) + T 3 M (3) + . . . j + Tipi + T 2 i> 2 + T 3 ^ 

= Uo + TAi + T 2 X 2 + T 3 X 3 + 0(T 4 )^ U + Tfa + T 2 ^ 2 + T 3 ^ (47) 

gives the following set of equations 





- *o)ipo 


= 








(m(°) 


- *o)ipi 


= (Ai 


-M«)0o 






(m<°) 


- *o)l/>2 


= (Ai 


-M^)ipi- 


h(A 2 




( M (°) 


- Ao)^ 


= (Ai 


-M«)^2- 


h(A 2 


-A/ (2) )^i + (A 3 -M^)^ 



Because ilf' 1 ) = 0, the second equation tells us that Ai = 0. This in turn means 
that ipi is an eigenfunction of Using the freedom to add any multiple of 

ipo to -0 S , s > 0, we can arrange (ipo, ips) = for all s > |S]. In particular, we 
may choose ip\ = 0. Hence, contracting the last equation with ipo gives us 

A 3 (Vo|Vo) = <Vo|M(%o>. (49) 
Going back to our notation 

n — 1 n — 1 

E^=E M ^--^ R oo. (50) 

i=l i=o 

However, as discussed earlier, the acceleration of a null-geodesics at q, though 
sub-leading in 0{T 2 ) is not sub-leading to 0(T 3 ) and hence cannot be ignored. 
It is at present not clear to us how to evaluate this contribution to the volume 
due to the acceleration. Instead, we resort to the Gibbons-Solodukhin approach 
by calculating the volumes of small causal diamonds in FRW spacetimes for 
which 8qR and <9o-Roo do not vanish. 
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Based on the considerations above we take the 0(T 3 ) correction to the vol- 
ume to be of the form 

V (x(n)d R(0) + K(n)d o Roo(0^j , (51) 

where again x( n ) an d /c(n) are universal constants to be evaluated. A more 
covariant version of the order T 3 terms is T a T 2 d a R(0) and T a T b T c d a R bc (0) 
which reduces to the above when the proper time is aligned along the time axis. 

We now calculate x( n ) an d k(u) by finding the volume of a small diamond 
in n-dimcnsional FRW spacetimes. It suffices to consider a spatially flat class 
of FRW spacetimes 

ds 2 = -dt 2 +t 2a ^{dx 1 ) 2 , (52) 
i=i 

where we have taken the scale factor to be of the form a(t) = t a . The Alexan- 
drov interval we will evaluate is centered at some t = to, with p and q chosen 
appropriately. We will examine two special cases a = 1 and a = 1/2. 



Case I: a = 1. 



We will work in conformal time, r\ = hi(^) or t = toexp I? , where we have 
chosen the integration constant so that r/ = at t = to, and to > 0. 



n-l 



ds 2 = 5(? ? ) 2 -drf + Y,( dxl ) 2 ( 53 ) 



where a{rf) = toe v . For this spacetime 

(n-l)(n-2) 
R(t) = f 2 -, i?oo = 

d R(t) = - 2(n ~ 1 ^ n ~ 2 \ 5 OJ Roo = 0, (54) 

which allows us to calculate x{ n )- 

We now calculate the volume of an interval I[p, q] where in conformal co- 
ordinates p = (— y , 0, . . . , 0) and q = , 0, . . . , 0), so that the proper time T 
between p and q is: 

T = t q — t p = t (e^ - e-T) = 2t sinh( — ), (55) 

or 

" = 2s "' h "<!;> = £ (1 -5;;r +0 < T4 »- (56) 
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Since g a b is conformally flat the boundary of I[p, q] is determined by the equation 



V(t) = V+ + V- 

fdr) a{v)" I dr r n ~ 2 I dfl n -2 + I dr] a(rj) n / dr r™" 2 I dfl n „ 



V f+v 

n—2 I jo i / J„ ;^\n / A n— 1 



2A -lt-jco S Hn V )^- v y 1 (57) 



The integral 



2 

f f N 

I cosh(ra7) I — — r] 



n iK I N\ , ,, I n - 1 



(2k)l 12/ v ; V I 

fe=0 z=o v ' v 7 v 



(I 



2-J n( 1+ 4(n + l)(n + 2) +Q ^)' 



where we have used the relation 



^~tT+^j{ I ) = a{ n-l ) ' (59) 



(60) 



1=0 

Using the expansion 

we find the expression for the volume 

V - Vv^-T 24(n+l)(n + 2)^ + ° (T ' 

- v °( i - T2 24(, + r ) (, + 2) ^ +o ( T4 < 61 > 

Clearly, even though d Q R(t ) = -2(n - l)(n - 2)t -3 ^ 0, the 0(T 3 ) terms 
vanish, hence implying that x( n ) = 0- 

Case II: a = 1/2. 
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Conformal time in this case is given by 77 = 2(sfi — \/to), where again 77 = 

2 

The conformal time scale factor is 



at t = t , tn > 0, so that t = f ( 1 + tAt- 

y 2 vto 

therefore 0(77) = a/^o + §■ ^ or * ms spacetime 



R(t) = 
doR(t) = - 



(n - l)(n - 4) 
it 5 : 
(n-l)(n-4) 



(n-1) 



doRoo{t) — — 



At 2 
(n-1) 



(62) 



Since we have already shown that x( n ) = 0; we can use this to calculate k(ji). 
For p, g given as before, the proper time T from p to q is: 

T = t q -t p = y/fcN, => TV = (63) 

Thus, 
V = V+ + V- 

2 _ 2 f ' 



(7?/ 



io + ^ 



/*(v5 + 2) / 



df7 



n— 2 



(n-1) 



in - 



AT 
T 



n-1 



„-l[n/2] 

n - 1) ^ 1 J 

v > l=o k=0 



2k 



drj rj 



l+2k 



(64) 



Using the relation Eq (|59[) and expanding in orders of N, this simplifies to 
2A„_ 2 fN^/tEYfl N 2 (n~l) 



V 



(n-l)V 2 J \n 16i (n + l)(n + 2) 

,2 !) , ^lATi\ 



0(N 4 ' 



V 1 + f 



16tg(n + l)(n + 2) 



0(7V 4 ) 



(65) 



Again, we see that the 0(T 3 ) terms are zero and we conclude that not only 
x(n) = but also k(ti) = 0. This suggests that although the term doRoo 
appears non-trivially in the tangent space Jacobian J, there is a non-trivial 
effect from the acceleration of the null-geodesics which must serve to cancel this 
contribution!]]- 

x Note that the 0(T 2 ) terms in both cases are compatible with our expression for the 
volume Eq (I34D to that order. 
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5 Conclusions 



In this note we have shown explicitly how to calculate the volume up to 0{T 2 ) 
for a small causal diamond l\p, q] in arbitrary dimensions using Myrheim's ap- 
proach. Our main result is that in any dimension the boundary term of Myrheim 
[T] takes the form 

/Ai[p ig ]=Vo(Jji*oo(0)Y (66) 

The dimension independence of the factor multiplying Vo is intriguing. At 
its most mundane, this stems from the fact that to this order it is only the 
sum of the eigenvalues of first order perturbation theory that contribute. This 
sum is itself dimension independent since it involves the trace J^^Ti -RoiOi(O) = 
S™=i -^oiOi(O) = i?oo(0). Whether there is a more profound reason for this 
universality is however unclear. 

Importantly, our analysis allows us to generalise this result to other integrals 
over the region I[p, q], so that for any integrablc function $(2;^) in I[p,q] 

J (Tx $(x") = ( 1 + ^j-Roo(O)) J d"z $(^) + 0(T 3+m ) (67) 

I[p><?] Io[p,<j] 

where m is the order of the flat spacctimc integral. This has useful implications 
for calculations in causal set theory. Recently it has been used to find an 
expression for the discrete Ricci scalar in a causal set in terms of the abundance 
of "k-chains" , a construction which differs substantially from that previously 
obtained by Bcnincasa and Dowker [9]. 

We have also extended our analysis to examine the order 0(T 3 ) contribu- 
tion. Using a combination of our analysis and the approach of Gibbons and 
Solodukhin [2] we find from calculations of the spacetime volume in two types 
of FRW spacetimes that the 0(T 3 ) contribution to the volume in fact vanishes. 
Since the spacetime volume plays such a key role in the continuum approxima- 
tion of causal set theory, we conjecture that the first non-trivial higher derivative 
correction to the Einstcin-Hilbcrt action from causal set theory must contain 
at least second derivatives of the curvature (modulo boundary terms) , thus dis- 
tinguishing the causal set effective action from other higher derivative theories. 
Our conjecture has obvious limitations since it rests solely on the assumption 
that these higher order corrections obtain only from contributions to the volume. 
We have neglected entirely the effect of other geometric and topological contri- 
butions that might arise in a more subtle fashion in causal set theory, although 
these are at present hard to construe. To substantiate the conjecture requires a 
far better understanding than we have at present of the discrete-continuum cor- 
respondence, and in particular the precise manner in which locality is recovered 
from causal set quantum gravity. 
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